import numpy as np

# 雅可比法的迭代矩阵 M
M = np.array([[0, -2, 2],
              [-1, 0, -1],
              [-2 , -2, 0 ]])


# 计算谱半径（最大特征值的绝对值）
eigvals_M = np.linalg.eigvals(M)


spectral_radius_M = max(abs(eigvals_M))


print(f"矩阵谱半径: {spectral_radius_M}")

